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ABSTRACT 

Gamma-ray bursts (GRBs) are presumed to be powered by still unknown central engines for 
the timescales in the range 1ms ~ a few s. We propose that the ORB central engines would 
be a viable site for strong meson synchrotron emission if they were the compact astrophysical 
objects such as neutron stars or rotating black holes with extremely strong magnetic fields 
H ^ 10^^ — lO^^G and if protons or heavy nuclei were accelerated to ultra-rclativistic energies 
of order ^ 10^^ — lO^^eV. Wc show that the charged scalar mesons like tt^ and heavy vector 
mesons like p, which have several decay modes onto tt^, could be emitted with high intensity 
a thousand times larger than photons through strong couplings to ultra-rclativistic nucleons. 
These meson synchrotron emission processes eventually produce a burst of very high-energy 
cosmic neutrinos with lO^^eV < E^. These neutrinos are to be detected during the early time 
duration of short GRBs. 

Subject headings: cosmic rays — elementary particles — gamma rays: bursts — magnetic fields 
— stars: neutron 
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1. Introduction 

Accumulating observational evidence shows the cxsistence of astrophysical objects with extremely 
strong magnetic fields ^ lO^^G. Kouveliotou et al. (1998) have revealed that a galactic X-ray pulsar with 
an estimated magnetic field of ~ 8 x lO^^G causes recurrent bursts of soft 7-rays, which are called soft 
gamma repeaters. Gamma-ray bursts (GRBs) with short and more intense bursts of lOOkeV ~ IMeV 
photons still remain puzzling although it has been progressively clearer that they are likely to have a 
cosmological origin. Some of the theoretical models of GRBs (Usov 1992, Kluzniak & Ruderman 1998, and 
Pacsynski 1998) invoke compact objects for their central engines such as neutron stars or rotating black 
holes with extremely strong magnetic fields H 10^^'"^''G in order to produce ultra-relativistic energy flow 
of a huge Lorentz factor T 10^. Furthermore, neutron stars or black holes with strong magnetic fields 
have been proposed by many authors (Hillas 1984, and references therein) as an acceleration site of ultra 
high-enrgy cosmic rays (UHECR). Milgrom & Usov (1995) suggested possible association of the two highest 
energy UHECRs with strong GRBs in their error boxes. 

An ultra-relativistic nucleus gives rise to efficient meson emission, in analogy with the canonical photon 
synchrotron radiation, in such a strong magnetic field because it couples strongly to meson fields as well as 
to an electromagnetic field. Due to its large coupling constant the meson synchrotron emission is ^ 10^ 
times stronger than the usual photon synchrotron radiation. Ginzburg & Syrovatskii (1965a, b) calculated 
the intensity of 7r° emission by a proton in a given magnetic field. At that time, however, they could hardly 
find the astrophysical sites to which their formulae could be applied. 

In the present paper we propose that GRB central engines could be a viable site for strong meson 
emission. Waxman & Bahcall (1997) proposed that high energy neutrinos with ~ are produced 

by photomeson interactions on shock-accelerated protons in the relativistic fireball (Rees & Meszaros 
1992). Pacsynski & Xu (1994) suggested that charged pions, which are produced in pp collisions when the 
kinetic energy of the fireball is dissipated through internal collisions within the ejecta, produce a burst 
of lO^^e^ neutrinos. Our proposed process of meson production is different from these two which can 
operate without magnetic fields. The very rapid variability time scale, 0.1ms, of many GRBs implies that 
each sub-burst refiects the intrinsic primary energy release from the central engine (Sari & Piran 1997). 
Thus, the length scale of the central engine is estimated to be 10km. The BATSE detector on board 
of Compton Gamma-Ray Observatory found the shortest burst of a duration of 5ms with substructure on 
scale of 0.2ms (Bhat et al. 1992). Accumulated BATSE data (Fishman & Meegan 1995) confirmed that 
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more than 25% of total events are short (< 2s) bursts. Therefore, if the central engines of GRBs were the 
compact stellar objects like neutron stars or rotating black holes associated with strong magnetic fields, 
relativistic protons or heavy nuclei would trigger meson synchrotron emission whose decay products could 
provide several observational signals even before the hidden explosion energy is transported to the radiation 
in the later stage of a relativistic fireball. 

We extend the previous treatment of 7r° emission of Ginzburg & Syrovatskii (1965a, b) to several kinds 
of neutral and charged mesons, which couple to a nucleon through scalar or vector type interaction, in a 
manner somewhat different from theirs. If the produced mesons are 7r°s or heavier mesons which have main 
decay modes onto 7r°s, they produce high energy photons which are immediately converted to lower energy 
7S or pairs through 7 — > 7 + 7 or 7 ^ e+ + e~ (Erber 1966) for their large optical depth in strong 
magnetic fields. However, if they are the charged mesons like tt^s, decay modes are very different and result 
in more interesting consequences. Since the source emitters of charged nuclei are accelerated to high energy 
(Hillas 1984) to some extent, very high energy neutrinos are produced by 7r+ — > /i+ + 1/^ — > e+ + f e + + t'ju 
or its mirror conjugate process for the similar mechanism as proposed by Waxman & Bahcall (1997). These 
neutrinos could be a signature for strong meson emission from the GRB central engines and would be 
detected in the very early phase of the short GRBs. 

We formulate the meson synchrotron emission in the next section, and the calculated results are shown 
and discussed in §3. 



We follow the established semi-classical treatment (Peskin & Schroeder 1995, Itzykson & Zuber 1980) 
of synchrotron emission of the quantum fields interacting with classical source current in strong external 
fields. In our present study the meson field is second quantized, while the nucleons are not and obey 
classical motion. 

The energy spectrum of synchrotron emission by a proton was first derived by Ginzburg & 
Syrovatskii (1965a, b) and it is give by 



where g is the strong coupling constant, /he w 14 (Sakurai 1967), E.,^ is the energy of the emitted pion. 



2. 



Semi-Classical Treatment of Meson Emission 




(1) 



and 7p is the Lorentz factor, 7^ = Ep/nipC 



? , of the rotating proton in a given magnetic field, -f^i/s is the 
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modified Bessel function of order 1/3. The function y{x) is given by 

where mpC^ = 938MeV and m,rC^ = 135MeV are the rest masses of proton and 7r° meson, and the 

parameter x^ which characterizes the synchrotron emission, is determined by the proton energy and the 
strength of magnetic field as x = jy^Ip with Hq = = 1.5 x lO^^G. In tlie above equations, variable x 
is introduced for mathematical simplicity, x — ^ x Since the available energy of the pion satisfies 
m-jr < < Ep, the corresponding variable range of a; is < a; < 

Applying the same treatment to vector mesons, we extensively obtain the energy spectrum of p meson 
synchrotron emission 

= (3) 

where K^/^ is the modified Bessel function of order 5/3. The function y{x) is defined by Eq. (2) by replacing 
m^r with p meson mass nipC^ — 770MeV, and x = ^ x Note that one can easily get the expression for 
photon synchrotron radiation in the limit of nip ^ by replacing the strong coupling constant g with the 
electromagnetic coupling constant e. 

The total intensity of a scalar or vector meson as a function of x is obtained by integrating equation 
(1) or (3) over meson energies m^_p < E,^^p < Ep or equivalently over 7"^ < a; < It is useful to give 

an approximate formula of the total intensity in the limit of large or small X'- 

3 r r \ (4) 

^^^Xexp(-f : X «1 



V3 fi' 



and 



^^(5/3)¥^x^/^ ^ ^ ^ : X » 1 

3/2\ 



x^^X^/^exp 



X <1 



where we have made an approximation Ki,{r]) w - — for 77 <C 1 (x S> 1), or K^{ri) « y^exp(— 77) for 
3> 1 (x ^ !)• These approximations are in reasonable agreement with exact numerical integrals within ± 
3% for TT meson and ± 10% for p meson at x < 0-01 or 10^ < x- 

Let us make a short remark on our classical treatment. When the proton energy is very high or the 
external magnetic field is strong, i.e. x ^ 1; the quantum effects not only in the meson field but in the 
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source nucleon current may not be negligible. In the case of photon synchrotron radiation, quantum effects 
were carefully studied (Erber 1966), and semi-classical treatment was found to be a good approximation to 
the exact solution within a few percent. In our treatment, the prefactor in Eq. (5) ( ^16^ ^(5/3)^ « 0.583. 
Taking the limit of rUp — > and g = e, Eq. (5) is applied to photon synchrotron radiation. It is shown 
(Erber 1966) that this factor should be 0.5563 in fully quantum mechanical calculation. It therefore is 
expected to hold true for the hadron processes as well. 



3. Results and Discussions 

Figure la shows a comparison between the calculated spectra of scalar 7r° meson emission and 7 
synchrotron radiation. Since Usov (1992), Kluzniak & Ruderman (1998), and Pacsynski (1998) suggested 
strong magnetic fields of order H 10^^'"^''G, which are presumed to associate with neutron stars or black 
holes of the GRB central engines, we here take H = 1.5 x lO^^G. The observed Lorentz factor of the 
fireball is F 10^, which indicates that the energy of charged particles is at least ^ Although 
the acceleration mechanism in GRBs is still unknown, there are suggestions (Milgrom & Usov 1995) that 
the GRBs associate with UHECRs. Six events of UHECRs beyond the Greisen-Zatsepin-Kuz'min cutoff 
energy ~ lO^^eF (Hill & Schramm 1985) have been detected by AGASA group (Takeda et al. 1998). We 
therefore vary the proton energy from Ep = to lO^^eF. All calculated spectra cut off sharply at 

incident proton energy E^^^^ = Ep. As seen in this figure, the intensity of tt^ emission for Ep = 10^^, 10^* 
and lO^^eV, which correspond respectively to x ~ 0.1, 10 and 1000, is 10^ ~ 10 times stronger than that of 
7 radiation in high energy parts of the spectra. This reflects the different coupling constants, (?^/e^ ~ 10^. 
Very sharp declination of the tt'^ spectra at lower energy arises from the integral in Eq. (1) for finite pion 
mass because K^{r]) ~ ^/^exp {—rj) for rj ^ 1, which is a good approximation in this energy region where 
y{x) » 1. 

Figure lb shows comparison between the calculated spectra of p meson emission and 7 synchrotron 
radiation. In this figure both spectra look very similar to each other, except for the absolute intensity and 
the sharp declination of low energy spectra due to finite p meson mass. This is because the interactions 
between p meson and proton and between photon and proton are of vector type. For the proton energies 
above lO^^eV the intensity of p meson emission is roughly a thousand times larger than that of 7 
synchrotron radiation, reflecting again its stronger coupling constant. 

Figure 2 displays the calculated total intensities of synchrotron emission of 7r° and p mesons as a function 
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of X- These are the integrated spectra shown in Figs, la & lb over available meson energies. The total 
intensity of 7 synchrotron radiation is also shown in this figure. I^o exceeds /-y at 3 x 10~^ < x < 3 x 10*, 
and Ip exceeds /-y at x > 0.2. Both I^o and Ip decrease exponentially with decreasing x due to their finite 
masses. (See the asymptotic forms at x ^ 1 in Eqs. (4) and (5).) This sharp declination of Ip takes place 
at higher x than I^o because p meson mass is larger than pion mass. Ip resembles I^ at x > 1, except that 
the intensity is a thousand times different from each other, reflecting that both p meson and photon have 
the same vector type coupling to proton with different couphng constatnts, g^/e^ ~ 10'^. 

A nucleon strongly couples to tt^ as well as field. The interaction Hamiltonian is charge independent 
and -ffint = ig{-J^i}n'yhi}p4>v+ + '4'p75'4'p4'tt°) + c.c. (Sakurai 1967). The initial state tjjp in both p ^ n + 7r+ 
and p — > p + 7r° processes is identical. The difference comes from the final states: Proton and 7r+ couple to 
external magnetic field, but neutron and do not. (Neutron has too small magnetic moment.) When we 
describe these processes in quantum mechanics, we should use wave functions from the solution of Dirac 
and Klein-Gordon equations for the nucleon and the pion separately. However, final states are more or less 
the same for the same charge state, aside from the different masses. These might not change the reaction 
amplitiudes by many orders at ultra-relativistic energies where the rest mass is neglected. Note that the 
conjugate processes n p + Tr~ and n — > n + 7r° can also occur when a composite nucleus like ^^Fe orbits 
in the strong magnetic fields. 

Heavy mesons including p meson have several appreciable branching ratios for the decay onto tt^. Let 
us discuss what kinds of observational signals they may make. 

We assume that some fraction of 10^^ ~ lO^^er^s of the gravitational or magnetic field energy is 
released by some unknown mechanism operating at the GRB central engine during very short time duration, 
~ 0.1ms, of the first sub-burst. We also assume that an appreciable part of this energy is deposited into 
the rclativistic motion of the material leading to UHECRs. In a somewhat different context, Waxman & 
Bahcall (1997) proposed that the photomeson production in the ejecta of the fireball would make a burst of 
neutrinos. Although their proposed mechanism of meson production is completely different from 
ours, we can apply similar discussion on the physical consequence of tt^ decay. Extending our previous 
discussions of in Fig. 2 to tt^, we expect that a thousand times stronger intensity of high energy 
neutrinos than the photons can be produced universally at 0.1 < x from Tr~^ ^ p.'^ + Vp ^ e'^ + + Up, + Up, 
and its mirror conjugate process. Since neutrinos can escape from the ambient matter, these neutrinos 
could be a clear signature showing strong meson synchrotron emission near the central engines of GRBs 
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associated with extremely strong magnetic fields. 

The neutron emerging from p ^ n + 7r+ inherits almost all proton energy and can escape from the 
region of strong magnetic field. If there is not a dense shell surrounding the central engine, it travels 
~ 10^ X (7„/10^") pc before beta decay. This process may also produce a very high energy neutrino. The 
generic picture of GRBs (Mcszaros & Rees 1993, Piran 1999) suggets that a baryon mass ~ 10~^Mq is 
involved in a single explosion. If this is the case, such amount of baryon mass is huge enough to stop almost 
all neutrons before running through the ambient matter. 

If the central engines are neutron stars or black holes, the material ejected from these compact objects 
contains heavy nuclei such as oxygen and iron because these are the products from evolved massive stars. 
Therefore, meson emission from a heavy nucleus as well as from a proton is worth being considered. Let us 

consider a nucleus of total energy E^tot, mass number A, and charge Z, in a magnetic field of strength H. 
The energy of each nucleon is ii^ = Etot/A. When the strength of the effective magnetic field is Hcs = , 
the orbital trajectory of a proton is the same as the trajectory of the nucleus in the magnetic field H . The 
intensity of tt" emission by the nucleus should be the sum of each nucleonic contribution provided that the 
synchrotron emission is incoherent. Thus the total intensity is given by 

(i?tot ,H)^Ax &^ {E, Feff ) , (6) 

where I^\e, Hes) is the intensity of 7r° emission by the nucleon of energy E in magnetic field Hgs . Note 
that both proton and neutron emit n^. Figure 3 shows 7^o^^(i^tot, H) for ^^Fe as a function of total energy 
Etot with a fixed magnetic field oi H = 1.5 x lO^^G. The sharp declination of 7^o^^ takes place at higher 
energy than 7^o^ because each nucleon in ^^Fe has effectively smaller energy than a single proton. 

There is now more motivation to study the GRBs in association with UHECRs. It is highly desirable to 
proceed with a project like Orbiting Wide-angle Light-collectors (OWL), in order to detect ultra-relativistic 
neutrinos from GRBs for finding the true nature of the central engines. 

This work has been supported in part by the Grant-in-Aid for Scientific Research (10640236, 
10044103,11127220) of the Ministry of Education, Science, Sports and Culture of Japan and also by 
JSPS-NSF Grant of the Japan-U.S. Joint Research Project. We thank the Institute for Nuclear Theory at 
the University of Washington for its hospitality and the Department of Energy for partial support during 
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5. Figure Captions 

Figure 1: (a) Calculated energy spectra of scalar n'^ meson emission (solid curve) and photon 
synchrotron radiation (dashed curve) for various incident proton energies Ep with a fixed magnetic field of 
ff = 1.5 X lO^^G. Denoted numbers in the figure are the proton energies Ep = 10^^, 10^^^, 10^'^, 10^*, 10^°, 
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and from left to right, (b) The same as those in (a) for vector p meson emission. 

Figure 2: Calculated total intensities of the emission of scalar tt" meson (solid curve), vector p meson 
(long-dashed curve), and photon 7 (dashed curve) as a function of x = 'ff^lp- 

Figure 3: Calculated total intensities of scalar -k^ meson emission by the proton (dashed curve) and iron 
nucleus (solid curve) as a function of the total energy Etot with a fixed magnetic field of H = 1.5 x lO^^G. 
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